GL-TR-90-0214 


A  Comparison  of  Manual  Scaled  and  Predicted 
foE  and  foFl  Critical  Frequencies 


Robert  R.  Gamache  and  William  T.  Kersey 


University  of  Lowell 

Center  for  Atmospheric  Research 

450  Aiken  Street 

Lowell,  Massachusetts  01854 


July  1990 


Scientific  Report  No.  2 


Approved  for  public  release;  distribution  unlimited. 


GEOPHYSICS  LABORATORY 

AIR  FORCE  SYSTEMS  COMMAND 

UNITED  STATES  AIR  FORCE 

HANSCOM  AIR  FORCE  BASE,  MA  01731-5000 


"This  technical  report  has  been  reviewed  and  is  approved  for  publication" 


M 

Co 


RGEN  BUCHAU 
ontract  Manager 


FOR  THE  COMMANDER 


This  report  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center,  all  others  should  apply  to  the  National  Technical  Information 
Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing  list,  or  if 
the  addressee  is  no  longer  employed  by  your  organization,  please  notify  GL/IMA, 
Hanscom  AFB,  MA  01731.  This  will  assist  us  in  maintaining  a  current  mailing  iist. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  requires  that  it  be  returned. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

ULRF-459/ CAR 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No.  0704-0188 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

University  of  Lowell 


6c.  ADDRESS  [City,  State,  and  ZIP  Code) 

Center  for  Atmospheric  Research 

450  Aiken  Street 

Lowell,  Massachusetts  01854 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Geophysics  Laboratory 


8c.  ADDRESS  [City,  State,  and  ZIP  Cods) 
Hanscom  Air  Force  Base 
Massachusetts  01731-5000 


6b.  OFFICE  SYM80L 
(If  applicable) 


8b.  OFFICE  SYMBOL 
(If  applicable) 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

GL-TR-90-0214 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Geophysics  Laboratory 


7b.  ADDRESS  (C/ly,  State,  and  ZIP  Code) 
Hanscom  Air  Force  Base 
Massachusetts  01731-5000 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F19628-90-K-0029 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 

124417F 


PROJECT 

TASK 

NO. 

NO. 

4643 

08 

WORK  UNIT 
ACCESSION  NO. 


11.  TITLE  (Include  Security  Classification) 

A  Comparison  of  Manual  Scaled  and  Predicted  foE  and  foFl  Critical  Frequencies 


12.  PERSONAL  AUTHOR(S) 

Robert  R.  Garaache;  William  T.  Kersey 


13a.  TYPE  OF  REPORT 

Scientific  No.  2 


16.  Sl!DPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

July  1990  56 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

foE,  foFl  predicted  critical  frequencies  , 

Ionospheric  modelling 


19.  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  CCIR  and  Titheridge  foE  critical  frequency  prediction  routines  were  tested  by  comparison 
with  1875  manually  scaled  values.  The  foFl  critical  frequency  prediction  routine  of 
Millraan  et  al  was  tested  by  comparison  with  1005  manually  scaled  values.  Plots  and 
statistics  of  the  comparisons  are  presented  and  discussed.  From  the  results  recommendations 
are  made  to  help  improve  autoscaling. 


20.  OISTRI8UTION/ AVAILABILITY  OF  ABSTRACT 
QUNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  OTIC  USERS 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Jurgen  Buchau 


DD  Form  1473,  JUN  86  Previous  editions  are  obsolete. 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22b.  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

LIS 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


TABLE  OF  CONTENTS 


Page 


1.0 

INTRODUCTION 

1 

2.0 

DATA 

3 

3.0 

PREDICTION  ROUTINES 

4 

4.0 

CALCULATIONS 

7 

5.0 

DISCUSSION  AND  CONCLUSION 

12 

6.0 

OTHER  GEOGRAPHIC  LOCATIONS 

15 

7.0 

REFERENCES 

44 

|  Acoession  Pop  j\ 

NTIS 
Di’IC 
Unanr 
J  ust  i 

GRAJfcl 

TAB  Q 

lounoed  □ 

fixation _ 

Bv 

Distribution/ 

Availability  Cedes 

foist 

Avail  and/or 
Special 

1 

iii 


LIST  OF  FIGURES 


Figure  No.  Page 

la  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  October  20,  1987,  Millstone  Hill,  MA  18 

lb  Error  foE(manuai)-foE(predicted)  vs. 

time  of  day  for  October  20,  1987, 

Millstone  Hill,  MA  18 

lc  foFl  (measured  and  predicted)  vs. 

time  of  day  for  October  20,  1987, 

Millstone  Hill,  MA  18 

Id  Error  foFl(manual)-foFl (predicted)  vs. 

time  of  dav  for  October  20,  1987, 

Millstone  Hill,  MA  18 

2a  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  February  20,  1988,  Millstone  Hill,  MA  20 

2b  Error  foE(manual)-foE(predicted)  vs. 

time  of  day  for  February  20,  1988, 

Millstone  Kill,  MA  20 

2c  foFl  (measured  and  predicted)  vs. 

time  of  day  for  February  20,  1988, 

Millstone  Hill,  MA  20 

2d  Error  foFl(manual)-foFl (predicted)  vs. 

time  of  day  for  February  20,  1988, 

Millstone  Hill,  MA  20 


IV 


LIST  OF  FIGURES  (Continued) 


Figure  No. 
3a 

3b 

3c 

3d 

4a 

4b 

4c 

4d 


Page 


foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 

for  March  3,  1988,  Millstone  Hill,  MA  22 

Error  foE(manual)-foE(predicted)  vs. 
time  of  day  for  Marc.i  3,  1988, 

Millstone  Hill,  MA  22 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  March  3,  1988, 

Millstone  Hill,  MA  22 

Error  foFl(manual)-foFl (predicted)  vs. 
time  of  day  for  March  3,  1988, 

Millstone  Hill,  MA  22 

foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 

for  March  4,  1988,  Millstone  Hill,  MA  24 

Error  'foE(manual)-foE(predicted)  vs. 
time  of  day  for  March  4,  1988, 

Millstone  Hill,  MA  24 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  March  4,  1988, 

Millstone  Hill,  MA  24 

Error  foFl(manual)-foFl (predicted)  vs. 
time  of  day  for  March  4,  1988, 

Millstone  Hill,  MA  24 


v 


LIST  OF  FIGURES  (Continued) 


Figure  No. 
5a 


5b 


5c 


5d 


6a 


6b 


6c 


6d 


Page 


foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 

for  March  5,  1988,  Millstone  Hill,  MA  26 

Error  foE(manual)-foE(predicted)  vs. 
time  of  day  for  March  5,  1988, 

Millstone  Hill,  MA  26 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  March  5,  1988, 

Millstone  Hill,  MA  26 

Error  foFl(manual)-foFl  (predicted)  vs. 
time  of  day  for  March  5,  1988, 

Millstone  Hill,  MA  26 

foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 

for  March  16,1988,  Millstone  Hill,  MA  28 

Error  foE(manual)-foE(predicted)  vs. 
time  of  day  for  March  16,1988, 

Millstone  Hill,  MA  28 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  March  16,1988, 

Millstone  Hill,  MA  28 

Error  foFl(manual)-foFl (predicted)  vs. 
time  of  day  for  March  16.1988, 

Millstone  Hill,  MA  28 


LIST  OF  FIGURES  (Continued) 


Figure  No. 
7a 


7b 


7c 


7d 

8a 


8b 


8c 


8d 


Page 


foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 

for  March  17,  1988,  Millstone  Hill,  MA  30 

Error  foE(manual)-foE(predicted)  vs. 
time  of  day  for  March  17,  1988, 

Millstone  Hill,  MA  30 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  March  17,  1988, 

Millstone  Hill,  MA  30 

Error  foFl(manual)-foFl(predicted)  vs. 
time  of  day  for  March  17,  1988, 

Millstone  Hill,  MA  30 

foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 

for  March  23,  1988,  Millstone  Hill,  MA  32 

Error  foE(manual)-foE(predicted)  vs. 
time  of  dav  for  March  23,  1988, 

Millstone  Hill,  MA  32 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  March  23,  1988, 

Millstone  Hill,  MA  32 

Error  foFl(manual)-foFl (predicted)  vs. 
time  of  day  for  March  23,  1988, 

Millstone  Hill,  MA  32 


V 


Vll 


I 


LIST  OF  FIGURES  (Continued) 


Figure  No.  Page 

9a  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  March  24,  1988,  Millstone  Hill,  MA  34 

9b  Error  foE(manual)-foE(predicted)  vs. 

time  of  day  for  March  24,  1988, 

Millstone  Hill,  MA  34 

9c  foFl  (measured  and  predicted)  vs. 

time  of  day  for  March  24,  1988, 

Millstone  Hill,  MA  34 

9d  Error  foFl(manual)-foFl(predicted)  vs. 

time  of  day  for  March  24,  1988, 

Millstone  Hill,  MA  34 

1 0  a  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  March  25,  1988,  Millstone  Hill,  MA  36 

10b  Error  foE(manual)-foE(predicted)  vs. 

time  of  day  for  March  25,  1988, 

Millstone  Hill,  MA  36 

10c  foFl  (measured  and  predicted)  vs. 

time  of  day  for  March  25,  1988, 

Millstone  Hill,  MA  36 

lOd  Error  foFl(manual)-foFl(predicted)  vs. 

lime  of  day  for  March  25,  1988, 

Millstone  Hill,  MA  36 


viii 


LIST  OF  FIGURES  (Continued) 


Figure  No.  Page 

1 1  a  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  April  12,  1988,  Millstone  Hill,  MA  38 

1 1  b  Error  foE(manual)-foE(predicted)  vs. 

time  of  day  for  April  12,  1988, 

Millstone  Hill,  MA  38 

12a  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  July  13,  1988,  Millstone  Hill,  MA.  40 

12b  Error  foE(manual)-foE(predicted)  vs. 

time  of  day  for  July  13,  1988, 

Millstone  Hill,  MA  40 

12c  foFl  (measured  and  predicted)  vs. 

time  of  day  for  July  13,  1988, 

Millstone  Hill,  MA  40 

1  2d  Error  foFl(manual)-foFl (predicted)  vs. 

time  of  day  for  July  13,  1988, 

Millstone  Hill,  MA  40 

1 3  a  foE  (measured,  CCIR  predicted,  and 

Titheridge  predicted)  vs.  time  of  day 
for  September  22,  1987,  Argentia,  NF  42 

13b  Error  foE(manual)-foE(predicted)  vs. 

time  of  day  for  September  22,  1987, 

Argentia,  NF  42 


LIST  OF  FIGURES  (Continued) 


Figure  No. 
13c 


13d 


Page 

foFl  (measured  and  predicted)  vs. 
time  of  day  for  September  22,  1987, 

Argentia,  NF  42 

Error  foFl(manual)-foFl(predicted)  vs. 
time  of  day  for  September  22,  1987, 

Argentia,  NF  42 


LIST  OF  TABLES 


Table  No.  Page 

1  Database  of  Manual  Scaled  Ionograms  3 

2  Measured-Predicted  Frequencies  in  MHz, 

AAD  and  AD  7 

3  Total  Statistics  for  the  Ionograms  Studied  11 

4  Statistics  for  Days  282  and  284,  1989  for 

Five  Stations,  Yearly  SSN  (158.9)  and 

10.7  cm  Flux  (212.5  16 

5  Statistics  for  Two  Months  at  other  Geo¬ 
graphic  Locations,  Yearly  SSN  and  10.7 

cm  Flux  16 

6  Total  Statistics  for  the  Ionograms  Studied  17 


1.0  INTRODUCTION 


The  ARTIST  autoscaling  routines  use  a  predicted  foE  to 
determine  a  range  to  search  for  the  ionospheric  characteristic  foE. 
When  there  is  no  measurement  of  foE,  the  predicted  foE  value  is 
reported.  There  are  several  formulations  that  can  be  used  to 
estimate  foE.  In  the  ARTIST,  the  predicted  foE  is  the  CCIR  model 
described  in  the  CCIR  Supplement  Report  252-2. 1  We  have  also 
tested  a  foE  prediction  routine  using  a  method  developed  by  John 
Titheridge.2  In  this  report  the  results  of  both  methods  are  compared 
with  the  manually  scaled  foE  for  1875  ionograms. 

For  determination  of  foFl  in  ARTIST,  the  search  range  is 
defined  at  1.2  to  1.9  times  foE.  In  the  future  this  range  will  be 
determined  by  a  foFl  prediction  routine.  We  have  tested  the  foFl 
algorithm  of  Millman  et.  al.3  The  results  of  the  method  were 
compared  with  the  manually  scaled  foFl  values  for  1,005  ionograms. 
A  program  to  predict  foFl  was  also  made  available  to  us  by  L. 
McNamara.4  However,  it  requires  a  file  of  numerical  coefficients  that 
is  used  by  interpolation  routines  that  are  too  slow  and  it  would 
occupy  too  much  memory  for  use  in  ARTIST. 

In  all  comparisons  two  differences  were  calculated,  the 
first  is  the  average  absolute  difference  AAD  defined  by 


AAD  = 


N 

£  |  Scaled  -  Predicted  | 

i  =  1 _ 

N 


0) 


and  the  other  is  simply  the  average  difference,  AD,  for  which 
cancelation  of  the  error  can  occur.  The  AD  is  useful  for  indicating 
when  prediction  routines  are  continuously  over  or  underestimating 
the  characteristics  whereas  the  AAD  is  a  better  indication  of  the 
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quality  of  the  routine.  One  should  keep  in  mind  when  studying  the 
statistics  that  the  manual  scaling  using  ADEP  yields  characteristics  at 
a  resolution  of  0.1  MHz. 

The  prediction  routines  require  the  solar  zenith  angle  for 
the  location  and  time  of  day.  The  CCIR  foE  prediction  uses  the 
sunspot  number,  SSN,  and  that  of  Titheridge  uses  the  10.7  cm  solar 
flux  number,  d>.  The  test  calculations  considered  the  yearly  mean, 
the  monthly  mean,  and  the  daily  values  for  SSN  and  <X>  for  each  day 
considered.  This  gives  an  indication  of  the  dependence  of  the 
predicted  critical  frequencies  on  these  parameters. 
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2.0  DATA 


Thirteen  days  of  data  were  chosen  for  the  initial  study, 
twelve  are  from  Millstone  Hill,  MA  and  the  other  is  from  Argentia, 
NF.  The  days  used  in  the  study  and  the  number  of  foE  and  foFl 
values  are  given  in  Table  1.  The  database  covers  magnetically  quiet 
to  moderately  active  days. 

Table  1.  Database  of  Manual  Scaled  Ionograms. 


Station 

Date 

#  foE 
values 

#  foFl 
values 

Millstone  Hill 

October  20,  1987 

99 

35 

Massachusetts 

February  20,  1988 

123 

67 

March  3,  1988 

42 

25 

March  4,  1988 

43 

29 

March  5,  1988 

43 

30 

March  16,  1988 

106 

81 

March  17,  1988 

84 

35 

March  23,  1988 

44 

29 

March  24,  1988 

46 

34 

March  25,  1988 

46 

29 

April  12,  1988 

July  13,  1988 

174 

121 

Argentia,  NF 

Sept.  22,  1987 

31 

22 

The  data  were  manually  scaled  using  the  ADEP  system 
and  files  of  the  scaled  characteristics  foE,  foFl,  and  foF2  were 
produced.  This  produced  a  database  of  911  foE  values  and  537  foFl 
values  to  test  the  prediction  routines. 
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3.0  PREDICTION  ROUTINES 


The  ARTIST  predicted  foE  is  the  recommended  CCIR 
method  described  in  the  CCIR  Supplement  to  report  252-2. 1  The 
value  of  foE  is  given  by 

foE  =  [(l  +  0.0094[d>i2  -  66])  cosm  Xnoon  (A  +  BcosX)  d]0,25  (2) 

where  O12  is  the  twelve  month  smoothed  solar  flux.  An  empirical 
relationship  exist  between  the  twelve  month  mean  sunspot  number, 
Rl2,  and  the  flux,<I>i2: 

0>i2  =  63.7  +  0.728  Ri2  +  0.00089  R2n  (3) 

In  equation  (2)  Xnoon  is  the  solar  zenith  angle  at  local  noon  and  X  is 
the  geographic  latitude  (positive  for  north  of  the  equator).  For  IX I  < 
32°  the  constants  in  equation  (2)  are  m  =  -1.93  +  1.92  cosX;  A  =  23 
and  B  =  116.  And  for  IXI  £  32°  the  constants  are  m  =  0.11  -  0.49  cosX; 
A  =  92  and  B  =  35.  D  is  the  time-of-day  factor  given  by 

(a)  for  %'  £  73° 

D  =  cosPx',  where  x'  is  related  to  the  solar  zenith  angle  x-  For  IXI 
£  32°,  x'  =  X>  but  for  IXI  >  32°  x'  is  taken  to  be  the  value  of  x  at  a 
time  0.05  hours  earlier.  For  values  of  IX!  <  12°  p  =1.31  and  for 
IXI  >  12°  p  is  given  the  value  1.20; 

(b)  for  73°  <  x'  <90° 

D  =  cosP  (x*  -  8x')>  where  8x'  =  6.27  *  10'13*(x'  -  50)8  in  degrees 
with  8x'  and  p  is  as  in  (a)  above; 
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(c)  for  %'  >  90  0 


D  =  (0.077)P  exp[-1.68  (tl-t)]  from  midnight  to  dawn  and 

D  =  (0.077)P  exp[-1.01  (t-t2)]  from  sunset  to  midnight 

where  t  is  the  local  time  of  interest  in  hours,  tl  is  the  local  time 
at  dawn  (x'  =90°)  in  hours  and  t2  is  the  local  time  at  sunset  (x' 
=  90°)  in  hours,  p  has  the  same  meaning  as  in  (A)  above.  We 
found  that  the  ARTIST  foE  prediction  routine  was  discontinuous 
from  (b)  to  (c),  i.e.  when  %'  becomes  greater  than  90°. 

The  other  foE  prediction  routine  is  taken  from  L. 
McNamara4  and  is  based  on  the  method  of  Titheridge2  (J.E.T.).  The 
formula  for  foE  is 

foE  =  1.12  •  4>0,25  •  cos(x)0,30  (4) 

where  for  x  >  70°,  cos(x)  is  replaced  by  the  inverse  of  a  Chapman 
function  evaluated  at  the  peak  of  the  E  layer  with  a  scale  height 
taken  as  YmE/2  and  zenith  angle  x> 

Both  methods  are  straight  forward  and  require  little 
computation  time.  The  12  month  smoothed  sunspot  number  and  the 
10.7  cm  flux  were  taken  from  the  National  Geophysical  Data  Center 
tables.5  The  comparisons  were  also  made  using  the  daily  and  mean 
monthly  values  for  SSN  and  <J>.  This  allowed  us  to  observe  the 
dependence  of  the  routines  on  these  values  and  to  compare  the 
quality  of  the  routines  using  daily,  monthly,  and  yearly  input  data. 
This  is  important  since  yearly,  and  perhaps  monthly,  values  are  the 
only  practical  values  to  run  the  algorithms,  as  they  are  used  in 
ARTIST. 


The  prediction  of  foFl  is  from  the  work  of  Millman, 
Bowser,  and  Swanson3  (MBS)  and  is  a  slightly  modified  version  of 
that  by  Davies.6  The  critical  frequency  of  the  FI -layer  is  given  by 
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foFl  =  (4.3  +  0.01  SSN)  cosn  (xfi) 


(5) 


where  XFI  is  the  normalized  FI -layer  solar  zenith  angle  given  by 

Xfi  =  90°  (x/105.5°)  for  0°  <;  x  ^  1 05.5°  (6) 

The  angle,  105.5°,  corresponds  to  the  zenith  angle  of  the  sun,  %  (in 
degrees),  when  illuminating  the  FI  layer  at  an  altitude  of 
approximately  240  km.  It  is  derived  by  assuming  the  absence  of  a 
screening  altitude. 

The  exponent  n,  in  equation  (5),  is  solar  zenith  angle 
dependent  according  to  reference  7 

n  =  0.2  for  0°  <;  x  *  90°  (7) 

n  =  0.2  +  0.3  (x  -  90°)/15.5°  for  90°  <  x  *  105.5°  (8) 

Equation  (8)  shows  that,  for  the  condition  90°  <  x  ^  105.5°.  the  value 
of  n  varies  linearly  between  0.2  and  0.5.  When  the  calculated  foFl 
value  at  a  location  is  greater  than  or  equal  to  model  foF2  values 
generated  from  numerical  fits  to  global  data,  foFl  is  assumed  absent. 

Note,  the  CCIR  prediction  routine  for  foE  and  the  foFl 
prediction  of  Millman  et  al.  use  the  SSN  while  the  foE  prediction  of 
Titheridge  uses  the  solar  10.7  cm  flux.  While  there  is  an 

approximate  relationship  between  the  two  [see  equation  (3)],  we 
have  chosen  to  use  flux  values  from  the  tables  directly.  Thus  in  the 
following  work,  the  mention  of  the  SSN  will  be  in  reference  to  CCIR 

prediction  routine  for  foE  and  the  foFl  prediction  of  Millman  et  al. 

and  the  mention  of  the  flux  will  be  in  reference  to  the  foE  prediction 
model  of  Titheridge. 
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4,0  CALCULATIONS 


Data  files  of  the  measured  foE,  foFl,  and  foF2  were 
created  with  the  station  location  and  times  of  the  soundings.  These 
were  input  to  a  program  which,  for  the  given  times  and  locations, 
calculates  the  predicted  values  of  foE  and  foFl.  The  routines  use  the 
SSN  and  the  solar  10.7  cm  flux  which  are  available  as  daily,  monthly, 
and  yearly  averages.  Thus  for  each  measurement  three  comparisons 
are  possible,  i.e.  for  predictions  using  the  daily,  monthly,  and  yearly 
averages  for  SSN  and  solar  10.7  cm  flux.  In  Tables  2  and  3 
comparisons  of  measurement  with  the  CCIR  predicted  foE  and  the 
Millman  et.  al.3  predicted  foFl  correspond  to  the  reported  SSN  and 
the  comparisons  with  the  Titheridge  model  predicted  foE  correspond 
to  the  reported  solar  10.7  cm  flux. 

The  difference  between  the  measured  and  predicted 
characteristics  was  formed  and  written  to  additional  files  for 
statistics  and  plotting.  In  the  tables  and  plots  all  times  are  reported 
in  UT,  the  differences  are  reported  as  Measured  -  Predicted,  and  the 
tables  report  both  the  AAD  and  the  AD.  Table  2  gives  the 
comparisons  for  the  thirteen  days  studied.  The  total  statistics  from 
the  911  measured  foE  values  and  the  537  measured  foFl  values  are 
shown  in  Table  3. 


Table  2.  Measured-Predicted  Frequencies  in  MHz,  AAD  and  AD 


SSN 

_ 

flunMI 

HSi 

foFl 

MBS 

Millstone  Hill  1988  dav  51  NfEl 

=  123 

NfFn  =  67 

Daily  Values 

51.00 

106.50 

AAD 

.09 

EH 

.28 

AD 

-.06 

H9 

-.21 

Monthly  Values 

40.00 

102.40 

AAD 

.07 

.10 

.22 

AD 

-.02 

-.11 

Yearly  Values 

100.20 

141.10 

AAD 

.25 

.16 

AD 

-.25 

-.13 

-.66  1 

Table  2.  Measured-Predicted  Frequencies  in  MHz,  AAD  and  AD 
(Continued) 


ne  Hill  1 


Daily  Values  72.00  99.10 

Monthly  Values  76.20  113.80 
Yearly  Values  100.20  141.10 


am 


nil 


Daily  Values 

77.00 

101.90 

Monthly  Values 

76.20 

113.80 

Yearly  Values 

100.20 

141.10 

.21 

.12 

-.21 

.05 

.21 

.10 

-.21 

-.02 

.30 

.18 

Daily  Values  64.00  102.60 
Monthly  Values  76.20  113.80 
Yearly  Values  100.20  141.10 


Daily  Values 

74.00 

114.10 

Monthly  Values 

76.20 

113.80 

Yearly  Values 

100.20 

141.10 

AAD 

AD 

AAD 

AD 

AAD 

AD 


AD  |-.19 
.20 
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Table  2.  Measured-Predicted  Frequencies  in  MHz,  AAD  and  AD 
(Continued) 


foFl 


Daily  Values 

99.00 

Monthly  Values 

76.20 

Yearly  Values 

100.20 

1  Millstone  Hill  1988  dav  8 

Daily  Values 

74.00 

Monthly  Values 

76.20 

Yearly  Values 

100.20 

I  BTC 


law 


Daily  Values  83.00  123.00 
Monthly  Values  76.20  113.80 
Yearly  Values  100.20  141.10 


illstone  Hill  1988  dav  85  N(E1  = 


Daily  Values  92.00  128.50 
Monthly  Values  76.20  113.80 
Yearly  Values  100.20  141.10 


=  29 


fan 


■ 

.15 

-.13 

.35 

JlJ. 


Table  2.  Measured-Predicted  Frequencies  in  MHz,  AAD  and  AD 
(Continued) 


SSN 

H | 

in 

foFl 

MBS 

103t  N(E)  =  77  ' _ ' _ 

Daily  Values 

118.00 

130.60 

AAD 

.44 

.38 

-- 

AD 

.23 

.35 

— 

Monthly  Values 

88.00 

123.60 

AAD 

.41 

.38 

AD 

.29 

.37 

— 

Yearly  Values 

100.20 

141.10 

AAD 

.42 

.38 

— 

AD 

,27 

.33 

-- 

Millstone  Hill  15 

>88  dav  195  N(E1 

=  174 

NCFl)  = 

121 

Daily  Values 

103.00 

141.30 

AAD 

.28 

.55 

AD 

-.06 

-.55 

Monthly  Values 

112.60 

157.60 

AAD 

.31 

.64 

AD 

-.09 

-.04 

-.64 

Yearly  Values 

100.20 

141.10 

AAD 

.27 

.18 

.52 

BiuSHI 

-.52 

Millstone  Hill  1987  dav  293  N(E) 

=  __99 

N(F1)  = 

35 

Daily  Values 

79.00 

95.60 

AAD 

.25 

.07 

.61 

AD 

-.25 

.06 

-.61 

Monthly  Values 

60.60 

97.40 

AAD 

.18 

.07 

.44 

AD 

-.18 

.05 

-.44 

Yearly  Values 

29.20 

85.30 

AAD 

.07 

.13 

.19 

■fll 

.13 

-.15 

1  Areentia.  NF  1987  dav  265  NflEl 

=  31 

NfFn  = 

22 

Daily  Values 

23.00 

81.30 

.08 

.17 

.13 

AD 

-.06 

.17 

-.13 

Monthly  Values 

33.90 

87.00 

AAD 

.12 

.13 

.23 

AD 

-.11 

.13 

-.23 

Yearly  Values 

29.20 

85.30 

AAD 

.10 

.14 

.18 

Hail 

Him 

_J4 _ 

-.18 

t  The  major  contributions  for  this  data  set  are  pre  to  sunrise  foE  values. 
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Table  3.  Total  Statistics  for  the  Ionograms  Studied 


SSN  and  10.7  cm  flux 


Daily  Values 


mmm 

AAD 

.18 

AD 

-.13 

AAD 

.17 

AD 

-.12 

AAD 

.22 

AD 

-.17 

Monthly  Values 
Yearly  Values 


.04 

.14 

-.05 


5.0  DISCUSSION  AND  CONCLUSION 


In  Table  2,  the  dependence  of  the  prediction  routines  on 
the  solar  indices  (sun  spot  number  or  10.7  cm  solar  flux)  can  be  seen. 
It  is  important  to  note  that  the  difference  in  the  errors  in  the 
predicted  characteristics  one  obtains  from  using  the  daily,  monthly, 
and  yearly  values  of  the  solar  indices  is  small.  These  results  are 
encouraging  since  from  a  practical  standpoint  the  yearly  or  monthly 
values  are  generally  used  in  the  algorithms.  It  appears  that  no 
penalty  jrs  from  using  the  yearly  mean  sun  spot  numbers  or 
solar  fl  values  in  the  foE  prediction  routines. 

For  each  day  studied,  plots  of  the  measured  foE  and  the 
two  predicted  values,  the  error  in  foE,  the  measured  foFl  and  the 
predicted  foFl,  and  the  error  in  foFl  were  made.  Figures  la- Id 
corresponds  to  the  first  entry  of  Table  1,  October  20,  1987  at 
Millstone  Hill,  and  the  figures  continue  in  this  manner  to  Figures 
13a-13d  the  last  entry  of  Table  1,  September  22,  1987  at  Argcntia, 
NF.  In  many  respects  the  figures  provide  more  information  about 
the  comparison  than  the  statistics,  giving  the  comparison  as  a 
function  of  time  of  day. 

Figures  were  produced  using  the  yearly  solar  indices 
(sunspot  and  solar  flux  numbers)  for  the  calculations.  In  the  figures 
1  through  13,  those  labeled  a  are  plots  of  the  measured  foE,  and  the 
two  predicted  values.  In  these  plots,  the  measured  data  are 
presented  by  the  dots,  the  CCIR  predicted  foE  values  are  given  by  the 
solid  line,  and  the  Titheridge  routine  predicted  foE  values  are 
represented  by  the  dashed  line.  The  plots  labeled  b  are  the  error 
given  by  foE(manual)-foE(predicted)  vs.  time  of  day.  The  CCIR 
comparisons  are  the  solid  circles  and  the  Titheridge  comparisons  are 
the  x  symbols.  The  figures  labeled  c  are  the  foFl  (measured-dots 
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and  predicted-solid  line)  vs.  time  of  day.  The  last  figure  per  group, 
labeled  d,  are  the  errors  in  the  predicted  foFl  values  vs.  time. 

In  comparing  the  results  for  predicting  foE,  on  all  days 
considered,  except  September  22,  1987  at  Argentia,  NF  and  October 
20,  1987  at  Millstone  Hill,  the  prediction  based  on  Titheridge's 
method  (the  dashed  line)  is  in  better  agreement  with  the  measured 
results.  Both  data  sets  that  have  CCIR  results  in  better  agreement 
are  from  1987.  This  is  probably  a  fortuitous  combination  of  data  set 
and  the  values  used  for  the  sun  spot  number.  In  general,  the 
predictions  track  the  measurements  rather  smoothly  but  are  usually 
higher,  but  not  always.  The  prediction  routine  of  Titheridge  appears 
to  track  the  manually  scaled  results  more  closely  than  the  CCIR 
prediction.  The  results  of  the  statistics  for  the  two  methods  indicates 
this  as  well.  The  averaged  absolute  difference,  AAD,  is  always 
smaller  for  the  Titheridge  routine.  More  important  is  that  the 
Titheridge  routine  is  more  centered  on  the  measured  values  as 
shown  by  the  average  percent  difference.  The  CCIR  routine  appears 
to  be  consistently  higher  than  the  measured  results,  this  can  be  seen 
by  comparing  the  AAD  to  the  AD  in  Table  3.  The  plots  of  the 
differences  (Figures  labeled  b)  do  not  offer  much  insight  into  the 
prediction  routines  except  supporting  the  previous  conclusion,  that 
the  Titheridge  predictions  follow  the  data  better  than  the  CCIR 
predictions,  leading  to  smaller  errors.  The  observed  structure  in  the 
error  plots  is  simply  due  to  the  digital  nature  of  the  measured  data. 
Based  on  the  results  it  would  appear  to  be  wise  to  adopt  the 
Titheridge  prediction  of  foE  for  use  in  the  ARTIST  algorithms. 

The  foFl  prediction  routine  gave  results  that  are 
generally  high  (see  Table  3).  Given  the  resolution  of  the  scaling, 
there  does  not  appear  to  be  a  significant  penalty  for  using  monthly 
or  yearly  mean  SSN  values.  By  comparing  the  AAD  with  the  AD,  one 
finds  the  foFl  prediction  to  be  consistently  higher  than  the  measured 
value.  This  is  demonstrated  in  the  plots  (set  c).  The  plots  also  show 
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that  the  foFl  prediction  does  not  track  measurements  as  smoothly  as 
the  foE  predictions.  The  foFl  prediction  may  be  quite  useful  for 
defining  a  search  range  for  the  autoscaling  routines.  From  the  foFl 
error  plots  (labeled  d),  we  note  that  the  search  range  defined  by 
foFl  predicted  -  1  -MHZ  to  foFlpredicted  +  0.2  MHz  contains  all  the 
measured  values  in  this  study. 
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6.0  OTHER  GEOGRAPHICAL  LOCATIONS 


Before  concluding  this  study,  the  prediction  routines  were 
tested  for  other  geographical  locations  for  which  data  were  now 
readily  available.  The  station  locations  and  statistics  from  the  study 
are  given  in  Table  4.  The  data  consist  of  two  days  where  recordings 
were  made  simultaneously  for  five  stations  as  shown  in  Table  4.  In 
addition  to  this,  data  was  also  available  for  one  month  at  Wallops 

Island,  VA  for  which  there  were  many  FI  measurements  and  for  one 
month  at  Kirtland  Air  Force  Base,  NM.  These  are  shown  in  Table  5. 

The  results  do  not  change  the  conclusions  from  above,  the  FI 

prediction  routine  gave  slightly  better  statistics. 

Adding  the  results  from  these  other  locations  to  the 
database  of  section  4  brings  the  total  number  of  foE  predictions 
studied  to  1875  and  the  total  number  of  FI  predictions  considered  to 
1005.  In  Table  6,  the  total  statistics  for  all  of  the  ionograms 

considered  are  reported.  For  the  additional  data,  the  daily  values  of 
the  sunspot  number  and  the  10.7  cm  flux  were  not  tested.  Thus,  the 
table  compares  the  monthly  and  yearly  values.  Again  these  results 
support  the  conclusions  of  section  5,  i.e.  that  there  is  no  significant 
penalty  for  using  yearly  or  monthly  mean  SSN  and  10.7  cm  flux 
values,  and  the  Titheridge  foE  prediction  gives  better  agreement  than 
the  CCIR  method. 
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Table  4.  Statistics  for  Days  282  and  284,  1989  for  Five  Stations, 
Yearly  SSN  (158.9)  and  10.7  cm  Flux  (212.5). 


fo 

CCIR  J 

E 

_ J.E.T 

foFl 

MBS 

Argentia,  NF  1989,  N(E)  = 
AAD 

AD 

n 

MSB 

mm 

.69 

.28 

Goosebay,  Labrador  N(E)  =  78 

AAD  |  .09 

AD _ t  -.08 

N(F1)  =  C 
.10 
-.05^ 

Bermuda,  FL  N(E)  = 

AAD 

AD 

46 
.06 
. -00 

N(F1)  =  1 
.09 
.03 

.03 

-.03 

Millstone,  MA  N(E)  = 

AAD 

AD 

■ 

HB 

1 

.61 

-.61 

Wallops  Island,  VA  N(S)  = 
AAD 

AD 

How 

JJSJRj 

1 

Table  5.  Statistics  for  Two  Months  at  other  Geographic  Locations, 
Yearly  SSN  and  10.7  cm  Flux. 


SSN 

CCIR 

foE 

1  J.E.T 

foFl 

_ MBS _ 

Wallops  Island,  VA 

N(E)  = 

553 

N(F1)  =  459 

July,  1988 

100.20 

141.10 

AAD 

.15 

.10 

.27 

AD 

-.11 

.01 

-.24 

Kirtland  AFB,  NM 

N(E)  =  43 

N(F1)  =  0 

January,  1990 

158.90 

212.50 

AAD 

.12 

.12 

AD 

.04 

.04 
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Table  6.  Total  Statistics  for  the  Ionograms  Studied. 


Solar  Indices 

foFl 

MBS 

..__N(EL= 

=  1875 

Monthly  Values 

AAD 

.17 

.11 

.38 

AD 

-.13 

-.36 

Yearly  Values 

AAD 

.17 

.12 

.40 

AD 

-.12 

-.37 
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Figure  la.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  October  20,  1987, 
Millstone  Hill,  MA. 

Figure  lb.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

October  20,  1987,  Millstone  Hill,  MA. 

Figure  lc.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

October  20,  1987,  Millstone  Hill,  MA. 

Figure  Id.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 

October  20,  1987,  Millstone  Hill,  MA. 
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Figure  2a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  February  20,  1988, 
Millstone  Hill,  MA. 

Figure  2b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

February  20,  1988,  Millstone  Hill,  MA. 

Figure  2c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

February  20,  1988,  Millstone  Hill,  MA. 

Figure  2d.  Error  foFl(manual)-foFl(predicted)  vs.  time  of  day  for 

February  20,  1988,  Millstone  Hill,  MA. 
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Figure  3a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  3,  1988,  Millstone 
Hill,  MA. 

Figure  3b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  3,  1988,  Millstone  Hill,  MA. 

Figure  3c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  3,  1988,  Millstone  Hill,  MA. 

Figure  3d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 

March  3,  1988,  Millstone  Hill,  MA. 
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Figure  4a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  4,  1988,  Millstone 
Hill,  MA. 

Figure  4b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  4,  1988,  Millstone  Hill,  MA. 

Figure  4c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  4,  1988,  Millstone  Hill,  MA. 

Figure  4d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 

March  4,  1988,  Millstone  Hill,  MA. 


Figure  5a.  foE  (measured,  CCIR  predicted,  and  Titheridge 
predicted)  vs.  time  of  day  for  March  5,  1988,  Millstone 
Hill,  MA. 

Figure  5b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  5,  1988,  Millstone  Hill,  MA. 

Figure  5c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  5,  1988,  Millstone  Hill,  MA. 

Figure  5d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 
March  5,  1988,  Millstone  Hill,  MA. 
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Figure  6a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  16,1988,  Millstone 
Hill,  MA. 

Figure  6b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  16,1988,  Millstone  Hill,  MA. 

Figure  6c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  16,1988,  Millstone  Hill,  MA. 

Figure  6d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 
March  16,1988,  Millstone  Hill,  MA. 
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Figure  7a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  17,  1988,  Millstone 
Hill,  MA. 

Figure  7b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  17,  1988,  Millstone  Hill,  MA. 

Figure  7c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  17,  1988,  Millstone  Hill,  MA. 

Figure  7d.  Error  foFl(manual)-foFl(predicted)  vs.  time  of  day  for 

March  17,  1988,  Millstone  Hill,  MA. 
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Figure  8a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  23,  1988,  Millstone 
Hill,  MA. 

Figure  8b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  23,  1988,  Millstone  Hill,  MA. 

Figure  8c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  23,  1988,  Millstone  Hill,  MA. 

Figure  8d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 

March  23,  1988,  Millstone  Hill,  MA. 
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Figure  9a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  24,  1988, 

Millstone  Hill,  MA. 

Figure  9b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  24,  1988,  Millstone  Hill,  MA. 

Figure  9c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  24,  1988,  Millstone  Hill,  MA. 

Figure  9d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 
March  24,  1988,  Millstone  Hill,  MA. 
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Figure  10a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  March  25,  1988,  Millstone 
Hill,  MA. 

Figure  10b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

March  25,  1988,  Millstone  Hill,  MA. 

Figure  10c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

March  25,  1988,  Millstone  Hill,  MA. 

Figure  lOd.  Error  foFl(manual)-foFl(predicted)  vs.  time  of  day  for 

March  25,  1988,  Millstone  Hill,  MA. 
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Figure  11a.  foE  (measured,  CCIR  predicted,  and  Titheridge 
predicted)  vs.  time  of  day  for  April  12,  1988,  Millstone 
Hill,  MA. 

Figure  lib.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 
April  12,  1988,  Millstone  Hill,  MA. 
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Figure  12a.  foE  (measured,  CCIR  predicted,  and  Titheridge 
predicted)  vs.  time  of  day  for  July  13,  1988,  Millstone 
Hill,  MA. 

Figure  12b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 
July  13,  1988,  Millstone  Hill,  MA. 

Figure  12c.  foFl  (measured  and  predicted)  vs.  time  of  day  for  July 
13,  1988,  Millstone  Hill,  MA. 

Figure  12d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 
July  13,  1988,  Millstone  Hill,  MA. 


Figure  13a.  foE  (measured,  CCIR  predicted,  and  Titheridge 

predicted)  vs.  time  of  day  for  September  22,  1987, 
Argentia,  NF. 

Figure  13b.  Error  foE(manual)-foE(predicted)  vs.  time  of  day  for 

September  22,  1987,  Argentia,  NF. 

Figure  13c.  foFl  (measured  and  predicted)  vs.  time  of  day  for 

September  22,  1987,  Argentia,  NF. 

Figure  13d.  Error  foFl(manual)-foFl (predicted)  vs.  time  of  day  for 

September  22,  1987,  Argentia,  NF. 
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